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f — ' After a discovery of the Higgs boson the next question is what are its couplings. At the LHC 

there should be many observable channels which can be exploited to measure the relevant 
parameters in the Higgs sector. Using the SFitter framework we map these measurements 
onto the parameter space of a weak-scale effective theory with free Higgs boson couplings. 
Our analysis benefits from the parameter determination tools and the error treatment used in 
new-physics searches, to study individual parameters and their error bars as well as parameter 
Oh' correlations. A special focus we will put on recent analyses using jet substructure techniques. 

<D ' 

1 Introduction 

^H ■ 

Understanding electro- weak symmetry breaking is one of the main goals of the LHC. In the 
Sta ndar d Model (SM) this is done spontaneously, and achieved by an SU(2) doublet, the Higgs 
field L^, which obtains a vacuum expectation value (vev). Three degrees of freedom become 
the longitudinal modes of W and Z bosons, while one remains as a physical scalar, the Higgs 
*i~ boson. The gauge boson masses arise from the kinetic term of the Higgs field in the Lagrangian 

Q | (D jU $)'''(D M $), which leads to WWHH and ZZHH terms, and replacing the Higgs field by its 

vev then to the masses. As these and the electromagnetic coupling are measured, this allows 
to determine the vev v = 246 GeV before actually observing the Higgs boson. Fermion masses 
we obtain from terms yfH^f^. The Yukawa couplings yf are a priori free parameters in the 
Lagrangian, but can be traded for the known fermion masses rrif = yt ■ v. This also means that 
the coupling of the Higgs to fermions scales as the fermion masses. 

The mass of the Higgs boson is the only unknown parameter in the SM. Electroweak precision 
tests tell us tha t it should be relatively light, just above the limit of 114.4 GeV from direct 
sear ches at LEP'^H Therefore we can use the theoretically predicted coupling values and test 
^jglojol a g a i ns t future LHC measurements'-^. Thereby we assume that the discrete quantum 
numbers, like spin and CP structure, of the boson are known^ and identical to the SM. Still 
there are many possible models which can generate such deviations. Examples inclu de sim ple 
extensions like ad ding another Higgs doublet, which is e.g. required in supersymmetry ' ', or 
composite models'^, where the Higgs emerges as a pseudo-Goldstone boson from a strongly- 
interacting sector. 

A correct treatment of errors is a crucial part for LHC parameter studies. Statistical errors 
from event counting are of the Poisson type. Systematic errors are correlated between individual 
measurements, so we need to include the full correlation matrix. Theory errors are best described 
as box-shaped using the RFit scheme'^. Using the SFitter tool ^ we can construct a fully- 
dimensional log-likelihood map of the parameter space. If we want to ask more specific questions, 



J for the SFitter collaboration 



Table 1: Signatures inclu ded i n our analysis for a Higgs mass of 120 GeV. The Standard Model event numbers 
for 30 fb _1 include cuts' — . The factor after the background rates describes how many events are used to 
extrapolate into the signal region. The last two columns give the one-sigma experimental and theory error bars 



on the signal. Table taken from Ref. 
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this must be reduced to lower-dimensional ones. Bayesian marginalization and Frequentist profile 
likelihoods are the two techniques we use, where the choice depends on the specific question we 
want to ask. 

In this article we will first review the LHC measurement channels which enter into our 
analysis. A particular emphasis we will put on the recently developed method of using jet 
substructure techniques^. Then we will present and discuss our results for Standard Model 
data. We focus on a Higgs boson of mg = 120 GeV, which is the preferred region from electro- 
weak precision data. 

2 Measurements 



At the LHC there are four main production modes of the Higgs boson: gluon fusion, weak-boson- 
fusion. associated production with vector bosons and associated production with a top-quark 
pairLyJ They need to be combined with the corresponding decay channels. For a light Higgs 
boson, like one with 120 GeV as we consider here, the main decay mode is into a pair of bottom 
quarks. Also decays via off-shell W and Z pairs lead to observable channels. Decays into taus 
can only be combined with weak-boson fusion^^, as we need to reconstruct the invariant mass 
of the tau pair, which is only known in the collinear limitM^. This channel is one of the discovery 
modes for a light Higgs and also allows us to determine its mass with a precision of ~ 5 GeV. 
The decay into photons is loop-induced and therefore only has a small branching fraction of 
~ 0.26% for a 120 GeV Higgs. Nevertheless the flat background, which is well subtractable by 
a side-band analysis, and a good 77 mass resolution make this mode a discovery channel and 
allow us to measure the mass with a precision O(100 MeV). 

In total we obtain the channels described in Table [TJ The numerical values for Hig gs pro- 
duction we get from Ref.'^ and for the decay from a modified version of HDECAY^M We do 
not include any channel which will only be measured at a later stage like the second-generation 
fermions. Such channels typically determine one additional parameter and therefore do not feed 



Figure 1: Illustration of the subjet algorithm. Starting with a fat jet, we first undo the last stages of the jet 

clustering until we have identified the bottom jets. We then apply a filter to remove the underlying event while 

retaining hard radiation from the Higgs decay products. Figure taken from Ref. 1 — " 

back into the analysis here. This also includes the Higgs self-couplings, which are important to 
establish the n ature of electro-weak symmetry breaking, but are notoriously hard to determine 
experimentally! I I 



3 Subjet Techniques 



In this section we will discuss in more detail the recent development of jet substructure tech- 
niques. Its presentation is based on the original paper^^. We consider associated production 
with vector-bosons with Higgs decays into a bottom-quark pair. A standard analysis would be 
overwhelmed by the large QCD background. Therefore we only consider a regime where both 
bosons are back-to-back and have large transverse momenta. For transverse momenta of the 
Higgs larger than 200 GeV only about one in twenty events survives. This is compensated by 
several advantages. The boost leads to central decay products, which can be tagged more easily. 
Also on-shell top quarks cannot simulate this kind of behavior and are no longer a background. 
Finally, the channel with Z decaying into neutrinos becomes visible as large missing transverse 
energy. 

Such a boost also means that often the jet reconstruction algorithm will not be able to 
resolve both bottom quarks. Instead it will combine them into a single fat jet. For R ~ ~^F 
this happens in roughly 75% of the cases. Therefore we need to identify such a fat jet. The 
Cambridge/Aachen jet algorithm-^ 1 has thereby shown to give the best results for the following 
procedure. Starting with a high-p? jet j we first undo the last stage of clustering so that 
we obtain two subjets j'i,j2 with rrij 1 > m,- 2 . Next we check if there has been a significant 
mass drop max(mj 1 ,mj 2 ) < 0.67m. Additionally we test if the splitting is not too asymmetric 

, ,„ > 0.09. If either of the conditions is not fulfilled, we take ?i as j and 



minfpE, . ,pL . ) ,.. 

— T ' 3 \ T ' J2 - AR? 



y 

.... 

repeat the steps. 

Otherwise we check if both subjets have b tags and if not, reject the event. This candi- 
date Higgs event is then filtered. We resolve the structure with a finer R separation i?gi t = 
min(0.3, R b i/2) and take the three hardest subjets. These are typically the two b jets and the 
leading radiation. A graphical description of the algorithm is depicted in Fig. [TJ 

The ATLAS collaboration has performed a full study using these techniques'^ and they 
obtain a statistical significance of 3.7a for a 120 GeV Higgs boson and a luminosity of 30 fb — 1 . 
Including 15% systematic uncertainty a significance of 3. Oct can be achieved. 

This method has also been applied to associated Higgs production with top quarks and bot- 
tom quark decays'-^. Here the subjet technique can help reduce the combinatorial background, 
which degrades the standard analysis. For a luminosity of 100 fb _1 the statistical significance of 
this channel is 4.1 standard deviations. Additionally it has been shown that further combining 
different subjet techniques can enhance the statistical significance even more^-A 



4 Calculational Setup 

For our analysis we assume a generalization of the SM Higgs sector with arbitrary couplings. 
Any coupling to particle j present in the SM we modify according to 

9jjH^g^ H (l + A jjH ) (l) 

where the Ajjh are independent of each other. As a global sign change in the Higgs couplings 
is not observable, we can take gw\VH to be always positive, or Awwh > — 1- Additionally 
loop- induced couplings to the photon and gluon are relevant. Here also the modified tree- 
level couplings enter. Also we allow for further dimension-five operators from new physics 
in the Lagrangian. An example for such a term are the additional loop contributions from 
supersymmetric partners. Therefore these couplings are modified to 

9jjh^9^h(1 + A.I¥h + &jjh), (2) 

where fffjjy is the loop-induced coupling in the SM, A^M^ the contribution from modified tree- 
level couplings and AjjH the additional dimension- five part. We also include the masses of 
the Higgs boson and the top- and bottom-quark into our parameter set. The errors on their 
measurements are large enough that their influence on the coupling determination should be 
taken into account. 

Furthermore we need to specify our treatment of the total width. This could also receive 
further contributions 

r tot =r t s ^(l + A r ), A r >0. (3) 

A simultaneous scaling of all couplings and the total width where g 4 /T tot = const, leaves all 
rates the same and we therefore cannot distinguish this. We will hence fix the total width to 
the sum of observed partial widths. 

The treatment of errors is a crucial part of the analysis. Statistical errors on individual 
channels are from counting experiments and therefore of Poisson type. The systematic errors 
are taken as fully correlated between the channels. They are derived from large event samples 
and can therefore safely assumed to be Gaussian. To combine these two types of errors we have 
devised an approximate formula of summing the inverse of the log-likelihoods 

i_^^_ = y 1 ( 4) 

X 2 -21ogL ^-21ogL i v; 

For a sum of Gaussians this gives the correct result of adding the errors in quadrature, while we 
have checked numerically that in all other cases we have a very good agreement with the exact 
way of a mathematical convolution. For theory errors we use the RFit prescription^^. For a 
deviation smaller than the theory error, the partial log-likelihood of the measurement is zero, 
while outside of this region it falls off with the combined experimental error. For the assumed 
numerical values of the errors we refer to Ref. . 

5 Results 

5.1 Likelihood map 

In Fig. [2] we show distributions for the three couplings Ay/WH-, A-ttH and A&&jj. In the left 
top row we show profile likelihoods for an unsmeared data set with an integrated luminosity 
of 30 fb _1 . Additional contributions to the loop-induced couplings we set to zero here. For all 
couplings we see a peak at the correct solution of A = 0. For both ttH and bbH a second peak 
at A = — 2 appears, which corresponds to a flipped sign of the coupling. For the top-quark the 
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Figure 2: Distributions for different couplings. Starting from the top left set with profile likelihoods for an 
unsmeared data set with an integrated luminosity of 30 fb _1 and no effective couplings, we modify these in the 
other sets: Bayesian marginalization (top right), 300 fb _1 (middle left), including effective coupling (middle right), 



and smeared data set with effective couplings (bottom). Plots taken from Ref. 
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likelihood of this solution is strongly reduced. The interference with a P^-loop in the effective 
Higgs-photon coupling allows us to determine the correct sign. In principle the same effect would 
also be true for the bottom quark. Its contribution to the effective photon or gluon coupling, 
which would fix the sign with respect to ttH, is too small to see an effect. 

On the right-hand side we take the same input in the top row, but now use Bayesian 
marginalization. The interference effect, which breaks the sign degeneracy in ttH, is completely 
washed out by volume effects. The bbH coupling shows a peculiar effect. Large values of the 
coupling are more likely than the correct one. The branching ratio into bottom quarks is much 
less affected by changing the coupling than other ones, because this also increases the total width 
significantly and the effect partly cancels. On the other hand large values allow larger changes 
in the other couplings without relevant changes of the overall rates. Therefore we have more 
parameter space available, so this is a pure volume effect. To understand the correlations we 
see Bayesian probabilities are less useful here and therefore we will show only profile likelihoods 
from now on. 

On the left-hand side of the middle row of Figure [2] we show a ten- fold increase of luminosity 
to 300 fb _1 compared to the top row. The general features stay unchanged, while the errors 
go down significantly, because the precision is primarily statistically limited. On the right-hand 
side we allow for additional contributions to the ggH and ^^H couplings at 30 fb _1 . We see 
that both peaks for the top-quark coupling now have the same height. The effective photon 
coupling can no longer break the sign degeneracy of ttH. Any mismatch is compensated by 
appropriately dialing the new term. 

In the bottom row we move from the true data set to a smeared one. This we have obtained 
by randomly smearing each measurement according to its errors. We see that the overall behavior 
does not change significantly. The best-fitting points move slightly away from their true values. 
Also the peak structure gets broader, as different measurements try to pull the parameters into 
different directions. 



Table 2: Errors on the measurements form 10000 toy experiments. We quote absolute errors on the couplings for 
30 fb , where additional contributions to the ggH and 77// couplings are either forbidden or allowed. For the 



latter we also show errors on the ratio of the coupling to WWH. Table taken from Ref. 
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5.2 Error determination 

To determine the errors we perform 10000 toy experiments, where we have smeared each mea- 
surement around the true data point including all experimental and theory errors. The resulting 
distribution we then fit with a Gaussian and extract cr symm . As the errors are not necessarily 
symmetric, we also fit two-half Gaussians with the same maximum and same height at the max- 
imum, but different errors below (<7 neg ) and above (<T pos ) the maximum. The fit-region is in all 
cases the central part within one standard deviation. 

From the results in Table [2] we see that the mass measurements hardly play any role. Their 
distributions are completely symmetric Gaussians and the errors correspond to the input values. 
The WWH coupling is the most precisely determined one. Its value does not change once we 
allow effective couplings, which means that the indirect determination via the loop-induced 
g^-yH is not important for a precise determination. The situation is completely different for 
ttH . As main contribution to g gg H and sub- leading to g 77 # the effective couplings significantly 
increase the error. Both main measurements for ttH and bbH are linked with a production 
mode where the WWH coupling enters, namely vector-boson fusion and associated production 
with subjet techniques, respectively. Therefore effective couplings do not change the errors here. 
The error on ZZH shows a particular effect. It decreases once we include effective couplings. 
This is an effect of the correlations between measurements. The main determination channel is 
gluon-fusion production with Higgs decay into a pair of Z bosons. Without effective couplings 
the production side is linked to the top-quark coupling, which is constrained by many other 
channels. Including them, this connection is removed and the error on the measurement need 
not be compensated by changing ZZH alone. Instead both can be changed, and due to the 
positive correlation between the two couplings this reduces the error. 

One might expect that forming ratios of coupling constants could improve the precision, as 
always certain combinations appear in our rate measurements. Due to its relatively small error 
the WWH coupling serves well as base value. Hence we also define a deviation on the coupling 
ratio to the WWH coupling as 



9jjH 
9WWH 



( 9jjH \ 
\9wwhJ 



SM 



1 + A, 



jjH/WWH 



)■ 



(5) 



The results we show in the right column of Table [21 In particular bbH benefits from forming 
the ratio. This is due to the total width, which appears in all measurements. The bbH coupling 



yields the largest contribution for a light Higgs boson and this leads to strong correlations to all 
other couplings. 

6 Conclusions 

In this article we have studied the determination of Higgs couplings at the LHC. As model 
we have assumed the Standard Model with free Higgs couplings, so that we are independent 
from the exact realization of new physics, if any. Using SFitter as a tool we map a set of 
LHC measurements onto the multi-dimensional parameter space. We have taken effects from 
statistical, correlated systematic and box-shaped theory errors into account. 

We find that we can determine the couplings with a precision of 20 — 40%. The improved 
accuracy from the newly developed subjet techniques is thereby an important ingredient. It helps 
to determine the bottom-quark Higgs coupling, which influences all others via its contribution 
to the total width. 
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